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Abstract 

We introduce and investigate the wellposedness of a model describing the self-propelled mo- 
tion of a small abstract swimmer in the 3-D incompressible fluid governed by the nonstationary 
Stokes equation, typically associated with the low Reynolds numbers. It is assumed that the 
swimmer's body consists of finitely many subsequently connected parts, identified with the fluid 
they occupy, linked by rotational and elastic Hooke's forces. In this paper we are attempting to 
extend the 2-D version of this model, introduced in [18j-[20j, to the 3-D case. Models like this 
are of interest in biological and engineering applications dealing with the study and design of 
propulsion systems in fluids. 
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1. Introduction and problem formulation. It seems that the first quantitative research in 
the area of swimming phenomenon was aimed at the biomechanics of specific biological species: Gray 
^0J(1932), Gray and Hancock [llj (1951), Taylor p6] (1951), [37J (1952), Wu [41j (1971), Lighthill 
[25| (1975), and others. These efforts resulted in the derivation of a number of mathematical models 
(linked the size of Reynolds number) for swimming motion in the whole R^- or i?^-spaces with the 
swimmer to be used as the reference frame, see, e.g., Childress [3J (1981) and the references therein. 
Such approach however requires some modification if one wants to track the actual position of 
swimmer in a fluid. 

A different modeling approach was proposed by Peskin in the computational mathematical 
biology (see Peskin [30j (1975), Fauci and Peskin [6J (1988), Fauci [7J (1993), Peskin and McQueen 
[31 J (1994) and the references therein), where a swimmer is modeled as an immaterial immersed 
boundary identified with the fluid, further discretized for computational purposes on some grid. 
In this case a fluid equation is to be complemented by a coupled infinite dimensional differential 
equation for the aforementioned "immersed boundary". 

In this paper we intend to deal with the swimming phenomenon in the framework of non- 
stationary PDE's along the immersed body approach summarized in Khapalov [20j (2005-2010). 
Namely, in [18j (2005), inspired by the ideas of the above-cited Peskin's method, we introduced a 
2-D model for ^^smair flexible swimmers assuming that their bodies are identified with the fluid 
occupying their shapes. This approach views such a swimmer as the already discretized aforemen- 
tioned immersed boundary supported on the respective grid cells, see, e.g.. Fig. 1 below. Our model 
offered two novel features: (1) it was set in a bounded domain with (2) governing equations to be 
a fluid equation coupled with a system of ODE^s describing the spatial position of swimmer within 
the space domain. We established the wellposedness of this model up to the contact of a swimmer 
at hand either with the boundary of space domain or with itself. The need of such type of models 
was motivated by the intention to investigate controllability properties of swimming phenomenon 
(see [20J). 

Our goal in this paper is to investigate the wellposedness of a 3-D version of this model. 
^This work was supported in part by NSF Grant DMS-1007981. 
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Figure 1: 4-parallelepipeds swimmer with all elastic forces and rotation forces about Z2{t) only. 

Further remarks on bibliography. It should be noted that, the classical mathematical issue 
of wellposedness of a swimming model as a system of PDE's for the first time was apparently 
addressed by Galdi [8| (1999) for a model of swimming micromotions in (with the swimmer 
serving as the reference frame). 

Another available approach to modeling of swimming motion (apparently initiated by the work 
Shapere and Wilczeck [3l] (1989)) exploits the idea that the swimmer's body shape transformations 
during the actual swimming process can be viewed as a set-valued map in time. The respective 
models describe swimmer's position via such maps, see fT3] (1981), [33j (2008), [4J (2011) and the 
references therein. Some models treat these maps as a priori prescribed, in which case the crux 
of the problem is to identify which maps are admissible, i.e., compatible with the principle of self- 
propulsion of swimming locomotion. In the case when the aforementioned motion map is not a 
priori prescribed (i.e., it will be defined at each moment of time by swimmer's internal forces and 
the interaction of its body with the resisting surrounding medium), the model will have to include 
extra equations, see, e.g., [40| in the framework of the immersed boundary method and the references 
therein. 

More recently, a number of significant efforts, both theoretical and experimental, were made to 
study models of possible bio-mimetic mechanical devises which employ the change of their geometry, 
inflicted by internal forces, as the means for self-propulsion, see, e.g., S. Hirose [16] (1993), Mason 
and Burdick [26] (2000); Mclsaac and Ostrowski [27J (2000); Martinez and J. Cortes [28j (2001); 
Trintafyllou et al. [39j (2000); Morgansen et al. [29j (2001); Fakuda et al. [5j (2002); Quo et al. 
[12] (2002); Hawthornee et al. [14j (2004), and the references therein. It was also recognized that 
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sophistication and complexity of design of bio-mimetic robots give rise to control-theoretic methods, 
see, e.g., Koiller et al. |[23| (1996); Mclsaac and Ostrowski |27j (2000); Martinez and Cortes [28] 
(2001); Trintafyllou et al. [39j (2000); San Martin et al. [32] (2007), Alouges et al. [Ij (2008), 
Sigalotti and Vivalda [35] (2009), and the references therein. It should be noted however that the 
above-cited results deal with control problems in the framework of ODE's only. 

A number of attempts were made along these lines to introduce various reduction techniques 
to convert swimming model equations into systems of ODE's, namely, by making use of applicable 
analytical considerations, empiric observations and experimental data, see, e.g., Becker et al 
(2003); Kanso et al. [H] (2005); San Martin et al. [32] (2007); Alogues et al. [Ij (2008), and the 
references therein. 

Problem formulation for 3-D swimming model. We consider the following model, con- 
sisting of two coupled systems of equations: one is a PDE system- for the fluid, governed by the 
nonstationary 3-D Stokes equation, and the other is an ODE system- for the position of the swim- 
ming object (or swimmer) in it: 

^ = uAy + F(z,v) - Vp in gT = f^x(0,T), y={y^,y2,ys) (1.1) 
div?/ = in Qt, y = in St = 50 x (0, T), y\t=o = yoi^^^ 

Zi{0) = Zio, z = l,...,n, n > 2, (1.2) 



,n, v{t) = {vi{t),...,Vn-2{t)) E R 



dt ~ mes{5,(0)} / 

S^{zi{t)) 

where for t e [0,r]: 

z(t) = (zUt), . . . , zJt)), zAt) G i?^ z = 1, . . 



F{z,v) = 5;][C.-i(x,t)fc. 



Zi{t) - Zi-l{t)\\ji3 -k-i) ... 
i-l ^-T7^ TTVn y^iyt) - Zi-i[t)) 



i=2 

+^i{x,t)ki-i- 



\\Zi{t) - Zi-l{t)\\R3 
Zi{t) - Zi-l{t)\\ji3 - k-l 



n-1 



i=2 
n-1 



\\Zi{t) - Zi-l(t)\\ji3 
C^-l(x,t) {Mz,_i{t) - Z,{t))) (t) - Cz+l(x,t) 



(z,_i(t)-z,(t))] 



Zi-i{t) - Zi{t)\\j^S . ... 

^ {Bi{Zi+l{t) - Zi{t))){t) 



Z^+l{t) - Z^{t)\\^ 



'^^i{x,t)Vi-i{t) 



i=2 



{Mz,_i{t) - Z,{t))){t) 



— — ^ {B,{z,^i{t) - Z,{t))){t) 



. (1.3) 



In the above, is a bounded domain in with boundary dVt of class C^, y{x,t) and p{x,t) are 
respectively the velocity and the pressure of the fluid at point x = (xi, X2, xs) G at time t, while 
z/ is a kinematic viscosity constant. Let us explain the terms in (1.1)-(1.3) in more detail. 

Swimmer: The swimmer in (1.1)-(1.3) is modeled as a collection of n bounded sets Si{zi{t))^i = 
1, . . . ,n of non-zero measure (such as balls, parallelepipeds, etc.), identifled with the fluid within 
the space which they occupy. These sets are assumed to be open bounded connected sets symmetric 
relative to the points Zi(t)^s which are their centers of mass. The sets Si(zi(t))^s are viewed as the 
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given sets 5^(0) 's ("0" stands of the origin) that have been shifted to the respective positions Zi{t)^s 
without changing their orientation in space. Respectively, 

!?!"ofcfL i = l,...,n. (1.4) 



0, if X e n\Si{zi{t)), 

Throughout the paper we assume that each Si{0) lies in a "small" neighborhood of the origin of 
given radius r > 0, while Si{a) denotes the set ^^(O) shifted to point a. Denote by 

5°= max {mes{5i(0)}}, 5o = min {mes{5i(0)}}. 

Forces: We assume that these sets are subsequently linked by forces described by the term F{z,v) 
No "actual" physical links between sets Si{zi{t)) are assumed (i.e., they are assumed to be negligible 
in terms of affect on the resisting surrounding fluid). The forces in (1.3) are internal^ relative to the 
swimmer - their sum is zero. We assume that a force applied to a set S{zi{t)) acts evenly upon all 
its points, and, as such, it creates an external force on the fluid surrounding S{zi{t)). 

The structural integrity of the swimming object is preserved by the elastic forces acting according 
to Hooke's Law. They act along the lines connecting the respective adjacent centers ^i(t)'s when 
the distances between any two adjacent points Zi-i{t) and Zi{t),i = 2,...,n deviate from the 
respective given values k-i > 0, z = 2, . . . , n as described in the first sum in (1.3). The parameters 
ki > 0, i = 1, . . . , n — 1 characterize the rigidity of the links Zi-i{t)zi{t)^ i = 2, . . . , n. The matching 
pairs of these forces between Zi-i(t) and Zi{t)^ and between Zi(t) and Zi^i{t) are shown on Fig. 1. 

The 2nd sum in (1.3) describes the rotation forces about any on the points ^^(t), z = 2, . . . , n — 
1 which make the adjacent points to rotate about it perpendicular to the lines connecting the 
respective Zi(t)'s. To satisfy the 3rd Newton's Law, these forces lie in the same plane along with 
the matching counter-force given in the 3rd sum in (1.3). Respectively, 

Ai = Ai{zi-i{t), Zi{t)), Zi^i{t)) : span {2:^_i(t) - Zi{t) , Zi^i{t) - Zi{t))}, 

Bi = Bi{zi-i{t),Zi{t)),Zi^i{t)) : sp^n {zi-i{t) - Zi(t), Zi^i{t) - Zi{t)} 

denote a non/mear mappings, defined at each moment of time by three vectors z^_i(t), ^i(t)), z^+i(t), 
such that for z = 2, . . . , n — 1: 



(z,_i(t) - z^{t)y[{A, {z,_i{t) - z,{t))){t)] = 0, {z^+iit) - z,{t)y[{B, {z,+i{t) - z,{t))){t)] = 0; 

II {Ai {zi-i{t) - Zi{t))){t) 11^3 = II Zi-i{t) - Zi{t) 11^3, 
II {Bi {zij,i{t) - Zi{t))){t) 11^3 = II Zi+i{t) - Zi{t) 11^3; 

• and the directions of vectors {Ai {zi-i{t) — Zi{t))){t) and {Bi {zi^i{t) — Zi{t))){t) are such that 
they correspond to either folding or unfolding motion of lines Zi-i{t)zi{t) and Zi{t)zi^i{t) 
relative to the point Zi{t). 

The magnitudes and directions of the rotation forces are determined by the given coefficients Vi{t)^ 
i = 1, ... ,n — 2. The choice of fractional coefficients at terms Ai{zi^i{t) — Zi{t)) in (1.3) ensures 
that the momentum of swimmer's internal forces is conserved at any t G (0, T) (see calculations in 
fT9] in the 2-D case). A matching pair of rotation forces, generated by Z2{t) for the adjacent points, 
is shown on Fig. 1. 
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Swimmer's motion. Dynamics of points Zi{t)^i{x,t),i = are determined by the 

average motion of the fluid within their respective supports Si{zi{t)ys as described in (1.2). 



Local and global approach to solutions of (1.1)-(1.3). Note that, when the adjacent 
points in the swimmer's body share the same position in space, the forcing term F in (1.3) and 
hence model (1.1)-(1.3) become undefined. While such situation mathematically seems possible, it 
does not have to happen. First of all, one can address the issue of local existence of solutions to 
(1.1)-(1.3) on some "small" time-interval (0,r), assuming that initially model (1.1)-(1.3) is well- 
defined in the above sense. This is the primary subject of this paper (see the next section). Then 
the question of global existence can be viewed as the issue of suitable selection of coefficients Vi^s 
with the purpose to ensure that the aforementioned ill-posed situation is avoided. 

In model (1.1)-(1.3) we chose the fluid governed by the nonstationary Stokes equation which, 
along with its stationary version, is a typical choice of fluid for micro-swimmers (the case of Low 
Reynolds numbers). The empiric reasoning behind this is that, due to the small size of swimmer, 
the inertia terms in the Navier-Stokes equation, containing the 1-st order derivatives in t and x, 
can be omitted, provided that the frequency parameter of the swimmer at hand is a quantity of 
order unity. However, it was noted that a microswimmer (e.g., a nano-size robot) may use a rather 
high frequency of motion, which may justify at least in some cases the need for the term yt in 
Stokes model equations. In general, it seems reasonable to suggest that the presence of this term 
(in a number of cases) can provide a better approximation of the Navier-Stokes equation than the 
lack of it. We also point out that in |6], [7], [3lj the full-size Navier-Stokes equation is used for 
micro-swimmers. It also seems that the methods we use for the nonstationary Stokes equation (as 
opposed to stationary Stokes equation), may serve as a natural step toward the swimming models, 
based on the Navier-Stokes equation. 

2. Main result: Local existence and uniqueness. Let J{ft) denote the set of inflnitely dif- 
ferentiable vector functions with values in which have compact support in Q and are divergence- 
free, i.e., dive/) = in fi. Denote by Jo{^) the closure of this set in the (L^(f]))^-norm and by G(fi) 
denote the orthogonal complement of Jo(^) in (^^(^))^ (see, e.g., [24], [38]). In J{Q) introduce the 
scalar product 

^33 

[01,02] = l^'^^ljx.hja^.dx, 01 (x) = (011,012,013), 02 (a:^) = (021,022,023). 

n 

Denote by H{Q) the Hilbert space which is the completion of J(fi) in the norm 



Il0i||if(^^) 



Everywhere below we will assume the following two assumptions: 

Assumption 2.1. For the given r > 0, defining the size of sets Si{0) in (1,4) ^ assume that 
/i_i > 2r, z = 2, Si{zi{jS)) d Vt, ||zi,o - %-,o||i^3 > 2r, i, j = 1, . . . , n, z 7^ j; (2.1) 

and the sets 5^(0), z = 1, . . . , n are such that 

[ dx = [ \Ux) - ^i{x - h)\dx < C\\h\\R3 yheBh,{0) (2.2) 

J (Si{o)us,{h))\{Si{o)nSi{h)) Jn 
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for some positive constants ho andC, where ^i{x) is the characteristic function of Si{0) and 5/^^(0) = 
{x III X ||^3< ho} C R^. 

Conditions (2.1) mean that at time t = 0, any two sets ^^(^^(O)) do not overlap, and that the 
swimmer hes in Jl. Condition (2.2) is a regularity assumption of Lipschitz type regarding the shift 
of the set 5'(0). It is satisfied, for instance, for balls and parallelepipeds. 

Assumption 2.2. Assume that within some (i?^)^ neighborhood G(z(0)) C (i?^)^ of the initial 
datum in (1.2) the mappings Ai and Bi are Lipschitz for a// i = 2, . . . , n — 1 in the following sense: 

II Ai{ai-i,ai, a^+i) (a^_i - ai) - Ai(hi-i, bi, 6^+i) {bi-i - bi)) ||^3 

< L{|| ai-i - bi-i) 11^3 + II - bi 11^3 + II a^+i - \\r^} , 
II Bi{ai-i,ai,ai^i) (a^+i - ai) - Bi{bi-i,bi,bi^i) (&^+i - bi) ||^3 

< L {II ai-i - bi-i) 11^3 + II a — bi ||^3 + || a^+i - 6^+1 ||^3} , 

for any a^±i, a^, bi±i^bi G G{z{0), where L > is a constant and Ai ^s and Bi ^s are defined as in the 
introduction by three respective vectors a^_i, a^, a^+i. 

Assumption 2.2 can be satisfied if, e.g., the points a^_i = 2;^_i(0),a^ = 2:^(0), a^+i = 2:^+1(0), i = 
2, . . . , n — 1 are not on the same line and the mappings Ai and Bi are selected by making use of 
Gramm-Schmidt orthogonalization procedure for vectors a^_i — ai and a^+i — a^. Alternatively, we 
can define A^'s and S^'s, making use of the cross-product: 

{Ai{zi-i{t) - Zi{t)){t) = ei{t) II Zi-i{t) - Zi{t)) 11^3, 

{Bi{zi^i{t) - Zi{t)){t) = e2{t) II Zi+i{t) - Zi{t)) 11^3, ei{t) = ,, ^if,, , i = 1,2, 

II \\r^ 

vi{t) = (z,-i(t) - z^{t)) X [(z^-i(t) - z,{t)) X (z,+i(t) - z,(t))], 

V2{t) = [(z^+l(t) - Z^(t))(z^+i(t) - Zi{t))] X (z,+i(t) - Zi{t)). 

Here is the main result of this paper. 

Theorem 2.1. Let yo G i?(f^); T > 0; ki > 0, i = \,...,n-\; Vi G L^(0,r), i = 
1, . . . ,n — 2; and 2:^(0) G 1^; i = 1, . . . , and let Assumptions 2.1 and 2.2 hold. Then there exists 
a T"" = r*(2;i(0), ... ,2:^(0), ||'?;i||i,oo(o^T); • • • ihn-2\\L^{o,T)M) ^ i^iT) such that system (1.1) - 
(1.3) admits a unique solution {y,p,z} on (0,T*), {y,Vp,z} G ^^(O, T*; 7^(1^)) x ^^(0, T*; G(l^)) x 
[C([0,r*];i?^)]-. Moreover, y G C([0, T*]; i/(f^)), yuy^,x,.Vx, e {L^Qt*))^ where ij = 1,2,3, 
and equations (1.1) and (1.2) are satisfied almost everywhere, while Assumptions 2.1 and 2.2 hold 
in [0,r*]. 

Remark 2.1. 

• The fact that conditions Assumptions 2.1 and 2.2 hold in [0,r*] implies that we are able to 
guarantee that within [0, T*] no parts of the swimmer's body will "collide", and simultaneously, 
that it stays strictly inside of fi. These conditions allow us to maintain the mathematical and 
physical wellposedness of model (1.1) - (1.3). 

• As it will follow from the proof below. Theorem 2.1 allows further extension of the solutions 
to (1.1) - (1.3) in time as long as Assumptions 2.1 and 2.2 continue to hold. This depends on 
the choice of parameters i;i(t), . . . , Vn-2{t)- 

Our plan to prove Theorem 2.1 is to proceed stepwise as follows: 
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• In Section 3 we discuss the existence and uniqueness of the solutions to the decoupled version 
of (1.2). 

• In Section 4 we will introduce three continuous mappings for the decoupled version of the 
system (1.1) - (1.3). 

• In Section 5 we will apply a fixed point argument to prove Theorem 2.1. 

In the proofs below we employ the methods introduced in [19j to investigate the wellposedness 
of the 2-D version of model (1.1)-(1.3), modifying and extending them to the ?>-D case. 

Without loss of generality, we will further assume that system (1.1)-(1.3) and all respective 
auxiliary systems below are considered on the time-intervals whose lengths are smaller than 1. 

3. Preliminary results. Introduce the following decoupled version of system (1.2): 

did) 1 C 

^ (Q (^w \ u(x,t)dx, wi(0) = Zi^o, z = l,...,n, (3.1) 

dt mes(S'^(0)) Js,{w,{t)) 

where u{x,t) is some given function. Denote w{t) = {wi{t), . . . ,Wn{t)). 

Lemma 3.1. Let T > and u G {L^{0,T; L'^{n))f be given. Then there is a T"" e (0,r) 
such that system (3.1) has a unique solution in C([0, T*]; i?^) satisfying Assumptions 2.1 and 2.2 
with w{t) in place of z{t), if they hold at time t = 0. 

Proof We will use the contraction principle to prove existence and uniqueness. Below the values 
of /lo, C are taken from (2.2). 

Select To to satisfy the following inequalities: 

['^ll'^ll(L2(g^))3 ^ ll'^ll(L2(0,T;L°°(17)))3 J 

< n^n min|7,;f'"""'°.^,'7'^'<'"»' T.ll. (3.2) 

Let for any given p G C([0, Tq];R^): 

Bho/2(p) -l^ze C([0,To];i?^) I \\z-p\\ci[o,To];R^) < y} ^ C([0,ro];i?^). 
For each i = 1, . . . , n, define a mapping Di : 5/^^/2 (^i,o) — > C{[0, Tq]; R^) by 

Di{wi(t)) = Zi,o H TFTFTW^ [ [ u{x,r)dxdr. 

mes(S'i(0)) Jo J s,{w,{r)) 

Then we can derive that: 



\\Di{wi{t))\\ji3 < Wzi^oWn^ + . ]\u\\fL2fQ^)Y V t G [0,To], z = l,...,n. (3.3a) 

Vmes(S'i(0)) 

Similarly, in view of (3.2): 

\\Di{wi{t)) - ^i,o||c([0,To];i?3) < ^=|==||w||(L2(Q^))3 < ^. (3.36) 
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Thus, Di maps Bf^^/2{^i,o) into itself for each i = l,...,n, where Zi^o is treated as a constant 
function. 

Let w^^\t) , w^'^\t) G S/^q/2(^z,o) and (^(x, 5) denote the characteristic function of a set 5 C R^. 
Then, making use of (2.2), we obtain: 

\\D,{wl'\t))-D,{wPm\ns = 

"t r ft 
u{x, T)dxdT — f 
Mi'Ht)) Jo JSi(wY>(r)) 



mes(5i(0)) 
1 

mes(5i(0)) 



/ / u{x, T)dxdT — / u{x, T)dxdT 

Jo JSiiwl'\t)) Jo JSi(wf\T)) 

1^ l^uix,T) {cix,Siiw\'\t)))-ax,SiiwP{T)))) dxdr 



fl3 



- mes(.gj(0)) ll^lk^^(0,r;L^ffi)))^ll^»-^^(^) ~ ^i^\i)\\c{[o,To]m^ [O'^o], i = l,...,n. (3.4) 

Therefore, after maximizing the left-hand side of (3.4) over [0,To], we conclude: 



llA(M^f^(t))-A(M^f^(t))||c([o,roi;i?3) < ^^^(^^^\M(L2(o,T;L^{n))r\\wi^^-wP (3.5) 
Hence, in view of (3.2), 

it fohows from (3.5) that Di is a contraction mapping on B 1^^/2(^1,0) ioi each z = 1, . . . , n. Therefore, 
there exist unique Wi{t) G (7([0, Tq]; i?^), i = 1, . . . ,n, such that Di{wi{t)) = i(;z(t), i.e., 

^2(^) = ^2,0 H TTTTTT^ / / ^(^^ T)dxdT, i = 1, . . . , n, (3.6) 

mes(S'^(0)) Jo JSi{wi{T)) 

which yields (3.1). 

On restrictions (2.3): Estimates (3.3a-b) imply that we may select a T* G (0,To) such that 
for any t G [0,r*], all Wi{tys will stay "close enough" to their initial values Zi^o^s to guarantee that 
Assumptions 2.1 and 2.2 holds for Wi{t)^ i = 1, . . . , n. This ends the proof of Lemma 3.1. 

4. Decoupled solution mappings. Let Bq{{)) denote a closed ball of radius q (its value will be 
selected in Section 5) with center at the origin in the Banach space 1/^(0, T; Jo(il)) H ^^(0, T; (i7^(Jl))^) 
endowed with the norm of 1/^(0, T; {H^{Q)f): 

B,{0) = {0 G L\0,T;Jomr\^\0,T; {H\n)f) \ Uh^o, T^iH^nm < q} , 

where H^{n) = {0 | 0, 0^^, 0^^^^. G L^ft), ij = 1,2,3}. 

Note that H'^{n) is continuously embedded into C(0), and thus L'^{0,T] {H'^{n)f) is continu- 
ously embedded into (L^(0, T; L^(ri)))^. This yields the estimate 



l(L2(0,T;L-(f7)))3 < JC||0||L2(O,T;(if2(Q))3) for SOmC K >0. (4.1) 

This implies that Lemma 3.1 holds for any u G Bq{0). 

4.1. Solution mapping for z^(t), i = 1, . . . , n. We now intend to show that the operator 
A:Bg{0) [C{[0,T];R')r, Au = w = (wi, . . . ,Wn), 
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where the Wi^s solve (3.1), is continuous and compact if T > is sufficiently small. 

Continuity. Let u^^^u^'^'^ G Bq{0) with Ti in place of T, where Ti > satisfies assumptions in 
the proof of Lemma 3.1 with Ti in place of T*. Define Au^^^ = w^^^ = {w[^\ . . . , w^"^) for j = 1, 2. 
To show A is continuous, we will evaluate 

term-by-term. To this end, similar to (3.4), we have the following estimate: 



77^— —V / / u^^\x,T)dxdT tt^t:^ / / u^'^\x,T)dxdT 

mes(^,(0)) 70 Js{wf\r)) mes(^,(0)) J S{wf\r)) 

= , ^ , II / / u^^\x,T)dxdT — f f u^'^\x,T)dxdT 

mes(5,(0)) 70 Js{w^\r)) Jo Js(wi'\r)) 



i?3 



+ 



< 



/ / u^^\x^T)dxdT — / / u^^\x^T)dxdT 

Jo Js{wI^\t)) Jo Js{wf\T)) 

^ I I {u^^\x,T)-U^^\x,T))£^{x,S{wf\T)))dxdT 

Jo Jn 



mes(5^(0)) 



i?3 



+ 



f I u('\x,rmx,S{wi'\r)))-ax.S{wi'\rm^^ 
Jo Jn 



< 



(1) - ,,.(2) 



< 



Vmes(5i(0)) 



' — u 



))3 + 



^^^^ll(L2(0,Ti;L-(n)))3|k)'^ - W^ric([0,Til;i?3 



(1) .,,(2), 



V'mes(S'o) 
Recall from (4.1) that 

So, select a T > as follows: 



mes(S'i(0))' 

«'"'ll(L2(0,ri;L-(f2)))3|kf^ -wP\\c([0,Ti];R3- (4-2) 



,(2)1 



mes(6'o) 



,(2)1 



(L2(o,Ti;L°°(a)))3 < 



< T < min 



/ mes(go) 
V CKq 



(4.3) 



Hence, replacing Ti in (4.2) with T satisfying (4.3) and maximizing the left-hand side of (4.2) over 
[0, r], we obtain: 



I (1) (2) II ^ 



(1) - ./(2) 



' — u 



ll(i2(Qj,))3 



^CKqVf (1) (2),, 



In view of (4.3), if \t) ^ on [0,r], then the above implies: 

( CKqVf \ (1) (2),, 



mes(S'o 



,(1) - „(2) 



9 



Thus, it follows that 



'^'^ - < v/rmes(go) _ ^(2)|| (4.4) 

\\C([0,T\,R ) - j^gg(^^) _ cKqVf^^ ^^^^^ ^ ^ 



Therefore, (4.3) and (4.4) imply that for every u^'^\ u^"^^ € Bq{Qi), 

IIA.(^) - A.(2)||[e([0,T];i.3)l. < - -(^^IU^(0,T;(^.(n))3). 

So, the operator A is continuous on Bq{Q) for sufficiently small T as in (4.3). 

Compactness. Furthermore, to show that A is compact, we will show that A maps any 
sequence in Bg(0) into a sequence in [C([0, T]; i?^)]^ which contains a convergent subsequence. To 
this end, consider any sequence {u^^^^^^^ in Bq{Qi). Using (3.6) with w\^^ and u^^^ in place of Wi 

and construct the sequence | , i = 1, . . . , n. 

Let us now show that is uniformly bounded and equicontinuous. Indeed, applying 

(4.1) to an estimate like (3.3a) and then maximizing over [0,r] yields: 

lkP^llc([o,T];i?3) < max {||^i,ob3} + . (4.5) 

*=iv..,ri ^mes(5o) 

To show equicontinuity, consider any t^t + h ^ [0? ^-g-? when h > 0. Then for i = 1, . . . , n: 



mes(5i(0)) 



Jt Js{wI'\t)) 



i?3 



^ ^ II (i)ii ^ • 1 

^mes(5i(0)) y^mes(5o) 

which implies the equicontinuity of , i = 1, . . . ,n (the case /i < is similar). Therefore, 

by Ascoli's Theorem, [Au^^^^^^ contains a convergent subsequence in [(7([0, T]; i?^)]^, i.e., A is 
compact on Bq{0). 

4.2. Solution mapping for decoupled Stokes equations. Now, consider the following 
decoupled Stokes initial boundary value problem: 

- uAy^ + Vp^ = /(x, t) in Qt, 

divy* = in Qt, = in St, y*|t=o = 2/0 ^ ^^(f^). (4.6) 

For any / G (L^(Qt))^, it is known that the boundary value problem (4.6) possesses a unique 
solution in L2(0, T; Jo(f^)) fl ^^(0, T; (i72(f^))3) ^ith the properties described in Theorem 2.1 
(see, e.g. [19], [27]). Moreover, (see, e.g. (7) on p. 79 and (50) on p. 65 in [19]), there is a positive 
constant L such that: 

lb*lli2(o,T;(i/2(Q))3) < L\\yo\\jj^^)+L / \\f{',T)\\%dT. (4.7) 

JQt 
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Thus it follows that, given yo, the operator 

B : iL\QT)f L\0, T; MQ)) f] L\0, T; iH\n)f), Bf = 

is continuous. 

4.3. The term ¥{z,v). Let T > be as in (4.3). Given u G 0^(0), let F^{w),w = {wi, . . . ,Wn) 
denote the value of F{z^ v) in (1.3), with Wi^s from (3.6) in z^'s, respectively. Consider the operator 

F: [C{[0,T];R^)r^{L\QT)f, Fw = F,{w). 

We will show that F is continuous, but first we evaluate ||F*(it;)||(^2(g^))2. To do that, we will use 
the standard algebraic transformation technique similar to that used for the 2-D swimming model 
in |19j . |20j . Therefore, we omit some of the details below. 
Let P{T) denote the upper bound in (4.5): 

P(r)= max {||^,o||^3} + ^l^=. (4.8) 

Then, using (4.5) with Wi in place of we can evaluate the first term in square brackets in the 
first line of (1.3) as follows: 

&(..r)t._ /ll"f -"■-''''ll«;-'---> («,W - 

\\Wi{t) - Wi-l{t)\\ji3 ^3 

< Co in, r, P(r), max {ki} , max {k}] , 

\ i=l,...,n—l i=l,...n—l J 

where and below Cq = Cq (n, r, P{T), max^^i^...^^_i {fc^} , max^^i^...^_i {U}) denotes a positive generic 
constant, continuous in r, fc^'s, /^'s, and P{T) (it can have a different concrete value for different 
expressions below). Along similar estimates for other terms in (1.3), we obtain: 

||F,H||(i2(Q^))3 < aV2^mes(0) (l + max {||t;i||ico(o,T)}') VnGi3,(0). (4.9) 

We will now show that F is a continuous operator. Let w^^^ ,w^^^ G [C([0, T]; i?^)]^. Consider 
(1.3) with wf^ in place of z^, z = 1, . . . , n, j = 1, 2. 

Making use of estimate like (3.3b), without loss of generality we can assume that T is small 
enough to ensure that Assumptions 2.1 and 2.2 holds (as stated in Lemma 3.1). Then, once again, 
similar to the %D case model in we can evaluate \F^{w^^^^ — F^{w^^^^\ji?> along the standard 
algebraic transformations and making use of Assumption 2.2, which will result in the following 
formula: 

WfM"-^) - FM^'^)\yL-{Q,)f < M(r)v/rmes(f2)||«;(i) - w^'^\\[ci[o,T]m]"^ (4-10) 

where 

M{T) ^ M{T,n,r,P{T), max {ki} , max {k})(l+ max {\\vj\\La.roT)}] , (4-11) 

i=l,...,n—l i=l,...n—l y j=l,...,n—2 J 

is defined similar to Co- Hence F is a continuous operator for T satisfying (4.3). 
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5. Proof of Theorem 2.1. Before proceeding to the proof of Theorem 2.1, we summarize 
the main resuhs of the previous section. In Section 4, we proved that for sufficiently smah T > 
satisfying (4.3), the operators 

A : B,iO) [Ci[0, T]; Au = w = (wi, Wn), 

F : [Ci[0,T];R^)r {L\QT)f , ¥w = F,{w), 

and 

B : {L\QT)f L^{0, T; J,(fi)) f| L^{0, T; {H\n)f), Bf = 
are ah continuous. A is also compact. As a result, 

BFA : Bq{0) L\0,T; J,(f])) fj L2(0, T; {H\n)f), BFAu = y, 

is continuous and compact. 

5.1. Existence: Fixed point argument. Select the value of q to be any positive number 
larger than V^||yo||i^^(Q)5 and choose T > as in (4.3) so that Lemma 3.1 holds. In view of (4.7) 
(with in place of /), we select T* G (0,min{T, 1}) small enough so that the continuous and 
compact operator BFA maps the closed ball Bq{0) into itself. 

Re-wrute (4.9) as follows: 

\mw)\\^LHQr))^<CiVf, 
where the positive constant Ci is independent of T. Select T* so that: 

o<r<.ni.. ^''y°' ,r,il. (5.1) 

Since T* satisfies (4.3), if we replace T with T* in (4.7), (4.9) - (4.11), we obtain: 

||BFA^t|||2(o,^*.(^2(Q))3) < L\\yo\\H{Q) + ^II^*('^)II(l2(q^*))3 

< L\\yo\\H(n) + LCfT'' < LWyoWj^^^^ + - L||2/o||^(q) = q^. 

Hence, BFA maps Bq{Q) into itself if (5.1) is satisfied. 

Thus, by Schauder's Fixed Point Theorem, BFA has a fixed point y which is a solution of the 
system (1.1) - (1.3), and which satisfies all of the requirements of Theorem 2.1. As usual, we may 
select Vp in L2(0,T*; G{9)) to complement the solution y e L^(0,r*; Jo{^)) in Theorem 2.1. This 
completes the proof of existence for Theorem 2.1. 

5.2. Uniqueness. To prove that the solution found in Section 5.1 is unique, we will argue by 
contradiction. Namely, suppose, e.g., that there are two different solutions 



and 



p) = (.f),...,4^)),y(^),p(^)} 



to (1.1)-(1.3), satisfying the properties described in Theorem 2.1 on some time interval [0, T], where 
T satisfies inequality (5.1). Without loss of generality, we assume the two solutions are different 
right from t = 0. 
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By (4.4), with in place of w^^\ i — 1, . . . , n, and y^-^^ in place of u^^\ both for j — 1, 2, we 
see that for any Tq G (0, T], and each z = 1, . . . , n: 

II (1) (2) II / v^Tomes(S'o) (2)11 /r 

Let us now evaluate Hy^-*^^ — y^^^ ||(l2(q^q))3. Note that (y^-*^) — y^'^^) satisfies the following Stokes 
initial- value problem: 

^^y"^"^ -y^'^^ = ^A(yW - + - F{z^^\v)) - V(pW in Qt,, 

div(y(i)-y(2)) = inQro, (^^'^ " y^'^) = in Eto, (y« - y('))|t=o = 0. 
According to (4.7) we have: 

U'^-y^^\\lLHQr,W r [ \\nz^'\v)-F{z^'\v)\\%dxdt. (5.3) 
In turn, similar to (4.10): 



for some N{T) is nonincreasing at T ^ 0^. Hence, combining (5.2) - (5.4) yields: 



\\F{z^^\v) - F{z^^\v)U. < N{T) J2 Ikf - #^llc([o,To];i?3). (5.4) 

i=i 



ll^yW-w(2)|L ^ niV(r)rov/Lmes(5o)mes(0) ^X^),, 

\\y y ll(L^(Q.o))3 < i^es(5(0)) - CKgVlb ^ ll(i^WTo))3- (5-5) 

Now, select Tq as follows: 

0<To<min| ^"^^^^^^ V ^ 4 • M 

\ 4niV(T) Vi^mes(5o)mes(f^) V ^CKq J j ^ ^ 

This choice of Tq implies that the following inequality holds: 

n7V(T)roA/^mes(S'o)mes(f^) + ^Ckq^/fo < ^mes(S'o). 
So, it follows from (5.5), (5.6) that: 

Therefore y^^^ = y^^^ on [0,To], and thus by (5.2), ^ = ^ for i = 1, . . . , n on [0,To]. Contradic- 
tion. This ends the proof of Theorem 2.1. 

6. Conclusion. In this paper we introduced a new hybrid model describing the locomotion of 
a "small" swimmer in the incompressible 3-D fluid. The model consists of two coupled systems of 
equations: one is a PDE system- for the fluid, governed by the nonstationary 3-D Stokes equation 
(typically associated with the low Reynolds numbers), and the other is an ODE system- for the posi- 
tion of swimmer in it. It is assumed that the swimmer's body consists of finitely many subsequently 
connected parts, identified with the fluid they occupy, which in turned are linked by the rotational 
and elastic Hooke's forces. We investigated the well-posedness of this model. Namely, in suitable 
function spaces we obtained the existence, uniqueness and regularity results for its solutions. These 
results are to be used for the derivation of a formula for asymptotically small motions of this type 
of 3-D swimmers, see (^^^ for the 2-D case). They can also be instrumental to study 

controllability properties of this swimming model [20j. Models like this are of interest in biological 
and engineering applications dealing with the study and design of self-propulsion systems in fluids. 
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